Generalization of geometric phase to completely positive maps 
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We generalize the notion of relative phase to completely positive maps with known unitary rep- 
resentation, based on interferometry. Parallel transport conditions that define the geometric phase 
for such maps are introduced. The interference effect is embodied in a set of interference patterns 
defined by flipping the environment state in one of the two paths. We show for the qubit that this 
structure gives rise to interesting additional information about the geometry of the evolution defined 
by the CP map. 
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Berry's |lj discovery of a geometric phase accompany- 
ing cyclic adiabatic evolution has triggered an immense 
interest in holonomy effects in quantum mechanics and 
has led to many generalizations. The restriction of adia- 
baticity was removed by Aharonov and Anandan S , who 
pointed out that the geometric phase is due to the curva- 
ture of the projective Hilbert space. It was extended to 
noncyclic evolution by Samuel and Bhandari || (see also 
Refs. ||]), based on Pancharatnam's f| work on inter- 
ference of classical light in distinct states of polarization. 
Another development of geometric phase was initiated by 
Uhlmann ^ who introduced this notion to mixed quan- 
tal states. More recently another mixed state geometric 
phase in the particular case of unitary evolution was dis- 
covered in the context of interferometry & . 

The geometric phase has shown to be useful in the 
context of quantum computing as a tool to achieve fault- 
tolerance H . For practical implementations of geometric 
quantum computing, it is important to understand the 
behavior of the geometric phase in the presence of deco- 
herence. For this, we generalize in this Letter the idea in 
to completely positive (CP) maps, i.e. we define the 
relative (Pancharatnam) phase and introduce a notion of 
parallel transport with concomitant geometric phase for 
such maps. These generalized concepts reduces to that 
of § in the case of unitary evolutions. 

Let us first consider a Mach-Zehnder interferometer 
with a variable relative U(l) phase x m one of the inter- 
ferometer beams (the reference beam) and assume that 
the interfering particles carry an additional internal de- 
gree of freedom, such as spin or polarization, in a pure 
state \k). The other beam (the target beam) is exposed 
to the unitary operator Ui, yielding the output interfer- 
ence pattern 7 oc 1 + ^cos(% — a), which is completely 



determined by the complex quantity 
ve ia = (k\Ui\k). 



(1) 



Thus, by varying x, the relative phase a and visibility 
v can be distinguished experimentally. We note that a 
is a shift in the maximum of the interference pattern, a 
fact that motivated Pancharatnam to define it as the 
relative phase between the internal states \k) and Ui\k) 
of the two beams. 

Pancharatnam's analysis was generalized in j?J to 
mixed states undergoing unitary evolution as follows. 
Assume that the incoming particle is in a mixed inter- 
nal state p = Ylk=i w k\k)(k\, where N is the dimension 
of the internal Hilbert space. Each pure component \k) 
of this mixture contributes an interference profile given 
by (fc|f7,-|A;) = i/ fe e lQfc weighted by the its probability Wk 
yielding 7 = J2k w kh oc 1 + J2k w k v k cos[x-a*]. Noting 
that Tr(Uip) = Wk(k\Ui\k), this can also be written as 
7 cx 1 + \Tr(Uip)\cos[x - argTr(L7 l p)]. The key result 
is that the interference fringes, produced by varying the 
phase x, is shifted by a — argTr(£7jp) and that this shift 
reduces to Pancharatnam's original prescription for pure 
states. These two facts are the central properties for a 
being a natural generalization of Pancharatnam's rela- 
tive phase to mixed states undergoing unitary evolution. 
Furthermore, it is clear that the quantity that extends 
that of Eq. (nl) to mixed states is 



ve 



(2) 



with visibility v = |Tr([/ip)|. 

Nonunitary evolution of a quantal state may be conve- 
niently modeled by appending an environment in a pure 
state that we designate |0 e ), i.e. 



Q = p®\0e)(0e 



(3) 
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FIG. 1: Interferometer for determining complete interference 
information. 



and letting the combined state evolve unitarily as g — > 
d = Ui e gU^ e with given [7j e . The evolved density ma- 
trix of the internal part is obtained by tracing over this 
environment yielding 



Tr e p' 



(4) 



where the Kraus operators are — (p e \Ui e \O e ) || in 
terms of an orthonormal basis {\p, e )}, p, = 0, . . . , K— 1 > 
N, of the A"— dimensional Hilbert space of the environ- 
ment. This map A is completely positive (CP), i.e. it 
takes density operators into density operators, and all 
trivial extensions I ® A likewise. Conversely, any CP 
map has a Kraus representation of the form Eq. (Q). 

Using Eq. (0) , the interference pattern for the incoming 
state in Eq. (||), evolved with {7j e in the target beam, is 
described by 



y e ia ° = Tr ie [U ie g] =Tr i [m p], 



(5) 



where we have used {O e \fx e ) = <5o M . The quantity «o 
is a natural definition of relative phase as it shifts the 
maximum of the interference pattern and reduces to the 
phase defined in (tJ for unitarily evolving mixed states. 

Since phase information has leaked from the system 
part, the interference information contained in Eq. (^|) is 
only partial. The remaining part may be uncovered by 
flipping the state of the environment associated with the 
reference beam to an orthogonal state |/x e ^ e ) (Fig. [I]) 
p2| . This transformation may be represented by the op- 
erator 



U = 




1 



<8> U it 



e lx 




U ® F . 



(6) 



where the first matrix in each term represents the spa- 
tial part and the operator F 0e ^ AIe flips |0 e ) to \(x e ^ e ). 
Fig. ^ shows the equivalent quantum network. The in- 
terference pattern is determined by 



= Tr^tWlOeKOelFo^J 

= r Tr le [U ie p\Q e ){p e \] = Tri[m M p] 



(7) 



for each fx — 1, ... ,K — 1. The set {i^e"*"}, p = 
0, . . . , K — 1, contains maximal information about the 




FIG. 2: Quantum network for interferometry of a state un- 
dergoing a CP map. 



interference effect associated with the CP map, by mea- 
suring on the system alone. For unitarily evolving mixed 
states one obtains i/„ = Sq^, due to orthogonality of the 
environmental states, and the surviving interference pat- 
tern Eq. (||) reduces to that of 0. 

The results above can be derived by considering purifi- 
cations. We may lift g in Eq. (||) to a purified state \^>) 
by attaching an ancilla according to 



|*) =^V^)|0e}|fca) 



(8) 



with {|fc a )} a basis in an auxiliary Hilbert space of dimen- 
sion at least as large as that of the internal Hilbert space. 
This state is mapped by the operators U — Ui e ® I a and 
F = U® F 0e ^ t i e <g> I a in the target and reference beam, 
respectively, i.e. 



|*tar) - V^[Uie\ki)\0 e )]\k a ), 
k 

|*rcf) - J2V^\k)[F 0e ^ e \0e}]\k a }, 



(9) 



Their inner-product becomes 



(*rcf|*tar> = ^ W k {h | <0 e | U ie |0 e ) | ki) 

k 

= y^ j Wk{h\m ll \ki) = Tr^m^p) (10) 
fe 

in agreement with Eqs. (||) and (f7|). 

To illustrate the above, let us consider the depolariza- 
tion channel JTo[ ] acting on a qubit in the initial state 
p = i(/ + r ■ a), where r = (x, y, z) is the Bloch vector 
with the length \r\ < 1, cr = (er x , a y , a z ) are the standard 
Pauli matrices, and / is the 2x2 unit matrix. We will 
model this with the Kraus operators 



(11) 



that map p — > p' = ^(1 + r' ■ cr). Here, mi, 7712, and 
7713, correspond to bit flip, both bit and phase flip, and 
phase flip, respectively, p is the probability that one 




3 



of these errors occurs and it determines the shrinking 
factor |r'|/|r| = (1 — 4p/3) of the Bloch vector. Exposing 
the depolarization channel to one of the interferometer 
beams, the interference pattern is determined by Eq. (|5|) 
as 

u e ia ° = T n (pm ) = VT^P- (12) 

This quantity is real and positive so that ao = 0, thus the 
channel only reduces the visibility by the factor y/1 — p. 
If the state of the environment is flipped to any of the 
other |/u e )'s the interference pattern Eq. ([?]) is determined 
by 

v^e**" = Tr i (pmJ = Mr li , (13) 

where r M = x,y,z for p, = 1, 2, 3, respectively. These are 
also real and positive so that = 0. The only effect is 
a nonvanishing visibility proportional to the probability 
amplitude \/p/3 associated with the corresponding er- 
ror. The absence of phase shifts can be understood from 
the fact that the depolarization channel only shrinks the 
length of the Bloch vector. 

Another illustration of the general formalism above is 
provided by the amplitude damping channel fiof , which 
models, e.g., the decay of an atom from one of its excited 
states to its ground state by emitting a photon. It may 
be described by the Kraus operators 

m a = ^[I + ^x) + g — ( /_CTz )' 

mi = ^-(ct x + icTy), (14) 

where p is the decay probability. When no photon has 
been emitted in the reference beam the interference pat- 
tern is determined by 

v e iao = Tt^pmo) = \ (l + VI - P + z[l - y/T=p\) , 

(15) 

which is real and positive so that ao = 0, and only the 
visibility v$ is affected anisotropically by the channel. 
On the other hand, if a photon has been emitted in the 
reference beam, we obtain 

^e Ml = Th(pmi) = ^y(* + *V), (16) 

being dependent only upon the x — y projection of the ini- 
tial Bloch vector. Here, the visibility decrease is propor- 
tional to the decay probability amplitude y/p and there is 
a shift in the interference oscillations determined by the 
angle of the x — y projection to the x axis. By checking 
whether or not there is a photon in the reference beam it 
should be possible experimentally to distinguish the two 
interference patterns determined by Eqs. ( |l5| ) and ( |l6| ) 
by post-selection. 



In the general case, it is convenient to make a polar 
decomposition of such that 

"fy = K u n> ( 17 ) 

where is Hermitian and positive, and is unitary. 
The action of each m M is uniquely defined up to N phase 
factors by the evolution of the system's density operator. 
This ambiguity must be associated with the correspond- 
ing unitarity as the Hermitian part is unique. 

Let us clarify this point for a qubit, where we may 
write — e r lB ^ e -^-' T anc j fo^ = a ^ -\- f,^ . a . it is possible 
to evaluate Eq. (Q) for p = \{I + r ■ er) yielding 

^e MfI = (a M + r ■ 6 M ) cos M + (r x e M ) • 6 M sin M 

-i(e M • bf, + a^r ■ e M ) sin M . (18) 

Now, with r = rn the density operator is unaffected 
by the change — > w p e~ l7f '™' <T as the additional part 
commutes with p. However, the interference pattern is 
determined by this new unitarity and therefore changes 
according to 

cos 9^ — > cos 9fj_ = cos 9fj_ cos 7 P — e p • n sin 9^ sin 7^ 
sin Ofj, — > e p sin 6> M — n cos M sin 7^ 

+e M sin 9^ cos 7 (I + e (l xn sin 6 M sin 7^(19) 

where u tl _e~ %1 ^ n ' rT — e - ,e f e M CT . 

The phase ambiguity can be removed by introducing 
a notion of parallel transport and concomitant geomet- 
ric phases for a continuous (time) parametrization of the 
CP map. For each interference pattern, the N parallel 
transport conditions read 

(k\ul (t)« M (t)|fc) = 0, k = l,...N, (20) 

where we have decomposed m At (t) = h^^u^it) — 
hfttyVpU^t) with w M = u M (0) completely specified by the 
channel, u M (0) = 1, and we have assumed h^(0) = S^q. 
These conditions naturally extend those in Eq. (13) of JtJ 
to the case of CP maps. They are sufficient and neces- 
sary to arrive at a unique notion of geometric phase in 
the context of single-particle interferometry. The set of 
these geometric phases for p = 0, . . . ,K—1 provides the 
complete geometric picture of the CP map in interferom- 
etry, given the unitary representation. 

Now, let us assume for simplicity that /i M and p diago- 
nalize in the same basis {|fc)} for all t > 0. In the cyclic 
case, where for each k we have u M : \k) — > e l ^ k \ k), /?£? be- 
ing the corresponding cyclic pure state geometric phase, 
the interference pattern is given by 

= ^v>k(k\h li \k)(k\v li \k)e i f , *\ (21) 
k 

where all (k\hn\k) > are real- valued as is Hermitian 
and positive. 



4 



Let us use Eq. ([21]) to compute the geometric phase 
for a qubit with r ^ in the depolarization channel 
with a unitary rotation added. Thus, we replace the 
Kraus operators m M in Eq. (|ll]) by m^u with it being an 
SU(2) operator that fulfills the parallel transport condi- 
tions Eq. (pp|). Here, is diagonal for all Kraus opera- 
tors, («i, V2, V3) = (a x , cr y ,a z ), and — u for all /1. For 
cyclic u, the interference patterns are determined by 



v e 

V 

vie 1 
v 3 e l 



/<I'd 



„i*3 



y/1 -p(cos(fi/2) + irsin(Q/2)) 





x/p/3{r cos(fi/2) + i sin(fi/2)) 



(22) 



with fl the solid angle enclosed by the loop on the Bloch 
sphere. The first interference pattern is precisely that ob- 
tained in |Q modified by a visibility factor ^/l — p. v\ and 
Vi vanish since the corresponding errors involve bit flips. 
Surprisingly, the pure phase flip in the last interference 
pattern introduces a nontrivial change in the appearance 
of r. This is due to the fact that this phase flip intro- 
duces an additional relative sign between the weights of 
the pure state interference patterns and is purely an ef- 
fect of the decoherence. 

We may also consider the case where u takes |0) — > 
|1) and vice versa. Any such u fulfilling the paral- 
lel transport conditions Eq. ( pp| ) may be written as 

e -i(n/2)[ca S <pa w +sm V <T y ]_ jj^ omy the errorg contain _ 

ing bit flips produce a nonvanishing interference effect 
that can be interpreted geometrically by noting that the 
SU(2) error operators iv\ and iv 2 are themselves taking 
|0) — * |1) and |1) — > |0) along geodesies intersecting the y 
and x axis, respectively. Thus, iv\u defines a closed loop 
on the Bloch sphere that encloses the solid angle 2ir — 2<p, 
yielding 



vie 



i'l'i 




(Q\iviu\0) 



1 - r 



\\iviu\\) 



— e ""^(cos <p + ir sin tp). 
3 



Similarly, iv 2 u defines the solid angle 37r — 2cp, so that 
the interference pattern is determined by 



v 2 e 



(<I>2 




iv 2 u\0) 



1 



(l|t«2tt|l) 



—e l7r (r cos tp + i sin ip) . 



(24) 



with known unitary representation, based on interferom- 
etry. We have further introduced parallel transport con- 
ditions that define the geometric phase for such maps. 
The interference effect is embodied in a set of interfer- 
ence patterns defined by flipping the environment state 
in one of the two particle beams. We have shown in the 
qubit case that this structure gives rise to interesting ad- 
ditional information about the geometry of the evolution 
defined by the CP map. We hope that this work will 
trigger new experiments on geometric phases for quantal 
systems exposed to environmental interactions. 
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Note added: After completing this work, it has come 
to our knowledge that Peixoto de Faria et al. |ll| have 
arrived at the interference pattern described by Eq. (|5|) 
in the context of measurement theory. 
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Again, there is a nontrivial change in the appearance of r 
in the last expression due to the fact that this error also rjji 
contains a phase flip. 

To summarize, we have provided a generalization of [12] 
the notion of relative phase to completely positive maps 
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This assumes we have full control over the environment 
and that it is, e.g., an extra degree of internal freedom. 



